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Abstract 

The main results of A. Zorich and I. Dynnikov about plane sections 
of periodic surfaces are extended to the PL case. As an application, the 
Stereographic Map of a truncated octahedron, extended to the whole R 3 
by periodicity, is analyzed numerically. 



1 Introduction 

The problem of the asymptotics of plane sections of smooth periodic surfaces, 
extracted from Physics literature by S.P. Novikov in 1982 Nov82 , turned out to 
be much richer then expected, leading ultimately to the association of a fractal 
on KP 2 to every element of a large class of triply periodic functions in K 3 (see 
Section 2). 

In order to visualize for the first time Stereographic Maps associated to triply 
periodic smooth functions and to get numerical confirmations of a Novikov con- 
jecture, claiming that the Hausdorff dimension of such fractals is stricly between 
1 and 2, we developed a C++ library and used it to investigate two smooth cases 
|Leo03| . Unfortunately, the running time of our numerical explorations grows 
way too much as soon as we sample with resolutions big enough to get a hint 
of the fractals, mainly because the number of polygons of the meshes approx- 
imating curved smooth surfaces (e.g. to retrieve its intersection with a plane) 
gets soon very big. 

This fact suggests that, from the numerical point of view, polyhedra are the 
best surfaces to study, at the very least for the obvious reason that the number 
of polygons needed to describe them at any resolution is constant. Moreover, 
such constant can be rather small even in non trivial cases, e.g. like in the case 
of the "extended truncated octahedron" , presented in this paper, which has just 
eight hexagonal faces. 
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Since the main results of the theory, due to A.V. Zorich |Zor84| and LA. 
Dynnikov |Dyn97| |Dyn99| , refer to the case of smooth surfaces only, in Section 
2 we will provide independent proofs for those theorems that make use of the 
smooth structure and will mention the main properties of the system. Then, in 
Section 3, we present the algorithm we implemented to explore numerically the 
problem and the results obtained in case of a the polyhedron obtained extending 
by periodicity the truncated octahedron. 

2 Fundamental objects and theorems 

2.1 Critical Points of height functions in polyhedra 

An analog of the Morse theory for height functions on polyhedra has been 
introduced by T. Banchoff in Ban67, Ban70;. Here we recall the concepts 
relevant for the present paper, slightly modified to cover the case of periodic 
polyhedra. 

Definition 1. By "embedded polyhedron" Mcl 3 we mean a countable collec- 
tion of cells K = {C r C K 3 } r= o,i,2j where 0-cells are points (vertices), 1-cells 
are closed connected segments and 2- cells are convex closed plane polygons, such 
that: 

PI the boundary of any cell is union of cells of lesser degree; 

P2 every cell having points in common with a higher degree cell is completely 
contained inside it; 

P3 K is locally finite, i. e. every vertex has a neighborhood in R 3 that intersects 
only finitely many cells; 

P4 for any point p € M , the union Star(p) of all cells containing that point 
is omeomorphic to an open disc. 

A triply periodic polyhedron is a polyhedron that is invariant with respect to a 
rank-3 discrete subgroup T ~ Z 3 of R 3 . Finally, by polyhedron M c T we 
mean the quotient M/T C M 3 /r ~ T 3 of a triply periodic polyhedron. 

Every triply periodic polyhedron M embedded in R 3 is the lift of a compact 
polyhedron M embedded in T 3 . Since we are going to study polyhedra's plane 
foliations, we are interested in the reciprocal relation between polyhedra and 
height functions (or, equivalently, constant 1-forms): 

Definition 2. A height function h(p) — h a p a is called "generic" for the poly- 
hedron M if no edge of M is perpendicular to the direction H = (h a ). Equiv- 
alently, a constant 1-form lo — h a dp a in R 3 (resp. T 3 ) is generic for M (resp. 
M) if no edge of M (resp. M) is contained in a single leaf 1 ofuj. 

1 Probenius theorem grants that the distribution u) = is integrable iff the 1-form u is 
closed. In this case, the leaves induced by u) = h a dp a in T 3 are the projections of the R 3 
planes perpendicular to H = (h a ); those induced on a polyhedron M. are the intersections of 
these leaves with the polyhedron. 
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Figure 1: Piece- wise linear Monkey Saddle: in the smooth however small 

perturbation would be enough to resolve this critical point in a pair of elemen- 
tary saddles; in the piece-wise linear case instead it is stable and, therefore, 
generic. 

Since height functions are not single-valued on T 3 (unless H is an integer 
direction, i.e. parallel to a lattice vector) while their differentials ui = dh are 
always well-defined in both K 3 and T 3 , we will refer mostly to 1-forms from now 
on. From the definition above it is clear that, like in the smooth case, the set 
of non-generic 1-forms has zero measure. 

Note that the foliation induced on T 3 (and therefore on M) by w does not 
change by multiplying the 1-form by a non-zero scalar, so from now on we will 
think of uj at the same time as a constant 1-form and as a point in MP 2 . 

Definition 3. The index of a point p 6 M (or, equivalently, [p] € M. ) with 
respect to a generic constant 1-form lu is the integer i(p,u>) = 1 — s/2, where 
s is the number of segments, having p as one of their extremes, in which the 
leaf of lu passing through p cuts Star(p). Ifi{p,uj) — 0, i.e. if the Star(p) is 
cut in exactly two components, the point is said "regular"; otherwise it is called 
"critical". 

The set of critical points for any generic constant 1-form is of course a subset 
of the set of vertices; exactly as in the smooth case, minima and maxima have 
index +1 and non-degenerate saddles have index —1. The main difference be- 
tween the smooth and the piece-wise linear (PL) case, for our purposes, is that 
saddles that are unstable in the former case, i.e. disappear for small perturba- 
tions of the 1-form direction, are stable in the latter and therefore cannot be 
disregarded; the simplest example is provided by the "Monkey saddle", which 
has index —2 (see Figure P). 

Nevertheless, an analog of the critical point theorem for generic 1-forms 
holds: 

Theorem 1. // M is a triply periodic polyhedron, invariant by the action of 
the rank-3 group T C Z , and to is generic for M , then Yl\p]eM *(Pj lj ) = x{M), 
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where [p] is the set of vertices I '-equivalent to p and the sum is hence extended 
to any set of inequivalent vertices. Equivalently, if we set to 1 the volume of a 
Dirichlet domain of T , then the average of the Euler characteristic of M con- 
verges to the Euler characteristic of M. : 

^pfj=lim J2 i(p,u)/Vol(B R ) = X (M) 

\\p\\<R 

Proof. Since Banchoff's proof Ba n70j of the critical point theorem for polyhedra 
is only based on local identities that are trivially true also in T 3 , that proof holds 
with no change for our case. 

The second part can be proved by considering that the genus g of the surface 
M. contained inside a cube of side R can be evaluated by reducing by homo- 
topy the surface to a graph and then evaluating the rank of the graph's first 
homology group. The result follows from the consideration that every compo- 
nent contained in a inner unitary cube contributes by g to the total genus of 
the component contained in the cube of radius R and their number grows with 
R 3 , while the cubes on the boundary provide a smaller contribute but can be 
disregarded in the limit for R — > oo since their number grows only as R 2 . □ 

2.2 Structure of foliations 

In the most general case, a constant 1-form lu induces on a triply periodic 
polyhedron M both open and closed leaves. 

Since being homotopic to zero is an open condition, leaves close enough to 
closed ones are also closed; maximal components of closed leaves are always 
enclosed between a pair of critical points of lu on M and form either cylinders 
(when both critical points are saddles) or discs (when one is a saddle and the 
other is a center) or spheres (when they are both centers). The last two cases are 
topologically trivial: a disc covered by closed leaves around a center is exactly a 
homotopy to a point of that component of M, and if M is a sphere then no open 
orbit can ever be induced on it by a closed 1-form. Hence, once the topologically 
trivial components are removed, what is left is a collection of cylinders Ci that 
separates a collection of subpolyhedra with boundary A/j filled by open leaves. 

Definition 4. A genus-k component of M is a PL submanifold with boundary 
M of M such that: 

1. dN is the finite disjoint union of plane parallel (topological) circles homo- 
topic to zero; 

2. the closed polyhedron AT, obtained by filling Af's holes with plane discs, 
has genus k. 

Definition 5. A Dynnikov decomposition Z of a polyhedron M. is a collection 
of subpolyhedra with boundary {Ci,J\fj} of Ai such that: 

1. every Ci is homotopic either to a closed cylinder or to a closed disc; 
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2. every Nj is a genus-kj component of M; 

3. i =/= j =>■ Mi l~l Mj = 0, while every other pair of distinct subpolyhedra of 
Z that is not disjoint shares single a boundary component; 

I M U,(UJ,A'.. 

The genus and rank of a Dynnikov decomposition Z are, respectively, the highest 
genus and rank of the MjS contained in Z; a genus- 1 rank-2 Dynnikov decom- 
position is called "Zorich decomposition". 

From the considerations above it is clear that every constant 1-form u in 
general position with respect to M, i.e. such that there are no saddle connec- 
tions, induces naturally a Dynnikov decomposition Z on M where all Cj are 
filled by the closed leaves and all TV} by the open ones. 

Under the same assumptions, to every focus it is associated a saddle that 
"cancels" it, namely there's a homotopy of the surface that gets rid of the pair 
saddle-focus without modifying the topology of the open orbits nearby; we call 
the saddle-type critical points that are left "topological saddles" , since it is them 
that contribute to the Euler characteristic of the surface. 

In the smooth generic case, the number of T-inequivalent topological saddles, 
i.e. the number of topological saddles in M, is of course exactly equal to 2c/— 2; in 
the PL case however, since we can have generic saddles with higher multiplicity, 
the Euler characteristic is only an upper bound for the number of topological 
saddles, as X (M) = 2-2g = £ feeZ k ■ #{\p]\i(w,p) = k}. 

2.3 The close-to-rational case 

Definition 6. By "irrationality degree" of a closed 1-form ft on a piecewise 
smooth manifold M we mean the number of rationally independent integrals ofCl 
over any base of the homology integer 1-cycles of M: irr(O) = diniQ spanqj J fl}, 
7 G Hi(M,Z); 1 -irrational forms are also called "rational". 

The structure of foliations induced on a periodic polyhedron by rational 
1-forms is much simpler than in the generic case since in this case all leaves 
are periodic; moreover, the appearance of a topological invariant enforces all 
1-forms close enough to rational to show the same behaviour. 

Rational 1-forms are rather special because foliations induced on M. by them 
are also induced by well-defined circle-valued functions on the polyhedron: in- 
deed, if {7i}i=i,2,3 is a base for i?i(T 3 ,R) and N an integer big enough such 
that N u) e Z, i = 1, 2, 3, and p a 6 M, then the function 

/ : M — ► S 1 

rP 

p i > exp(iN / ui) 

is well defined and differentiable, and its differential df p = iNcj p f(p) is pro- 
portional to uj; since / is never zero, the set df = on M coincides with the 
restriction to M of w = 0. 
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Figure 2: A (hypothetical) saddle point where four open leaves meet simultane- 
ously: the two branches of each of the inner open leaves cannot be adjacent to 
each other, so this picture cannot come from the section of a locally euclidean 
surface. 

This shows that all leaves induced by rational 1-forms on a polyhedron M. C 
T 3 are compact and therefore corresponding leaves on M C K 3 will be open or 
closed according to their homology class: leaves homotopic to zero in T 3 will 
remain so in any covering, all others will open up in T 3 's universal covering. 

Definition 7. From now on we will refer to closed leaves non homotopic to 
zero as "periodic" or, more generically, "open" leaves, so that by "closed leaf" 
we will implicitly mean a closed leaf homotopic to zero. 

Lemma 1. Be u) a rational 1-form in general position with respect to a polyhe- 
dron M: then no more than two open leaves can collide at any saddle point. 

Proof. The 1-form ui foliates T 3 in a 1-parameter family of embedded 2-tori, so 
that all open leaves at the same level of oj are parallel (i.e. they represent the 
same 1-cycle modulo sign); moreover the number of open leaves on every level 
is even, since they are all indivisible and their sum must be zero. 

It is easy to check, just by drawing pictures, that it is possible to have saddles 
of any index with closed leaves, and adding a single pair of open leaves does 
not change this situation; their presence though does not allow the presence of 
any other pair, since in any saddle point the two extremes of the same (critical) 
leaf must appear next to each other and this is of course impossible for all open 
leaves, apart for the two most external ones (e.g. see Figure [5}. □ 

Theorem 2. The Dynnikov decomposition induced on M by any constant ra- 
tional 1-form in general position is Zorich. 

Proof. Since all (non critical) leaves induced by w on X are circles, the critical 
points of lj on M. determine a subdivision of the polyhedron in the connected 
sum of a finite number of cobordisms. In the smooth case, the only non trivial 
cobordisms are elementary (i.e. pants-like) and Zorich Zor84 proved that in 
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the boundary of each pant there is at least one of the three boundary loops that 
is homotopic to zero (in T 3 ), from which the theorem follows easily. 

In the polyhedra case the cobordisms are not necessarily elementary, i.e. 
more than two leaves may collide at a saddle point, but lemma ^ grants us 
that even in this case no more than two loops in the cobordism may be open. 
This shows that the components of open leaves have genus 1 and therefore the 
decomposition induced on M. by to is Zorich. □ 

The presence of a genus- 1 rank-2 component of open leaves is a very strong 
condition and leads to the appearance of a topological quantity associated to 
the foliation |Dyn97| : 

Theorem 3. A Dynnikov decomposition Z of M. with rank 2 has genus 1 if 
and only if at least one of its Mj embedded with rank 2 has genus 1. In this case, 
all Mj are genus- 1 components of Mi and those embedded with rank 2 represent, 
modulo sign, the same indivisible non-zero 2-homology class in T . 

Definition 8. We call "soul" of a Zorich decomposition Z the "unsigned" non- 
trivial indivisible homology class I G PH2(T 3 , Z) common to all rank-2 compo- 
nents of Z . 

Corollary 1. Be lu a rational 1-form inducing on Mi a rank-2 Zorich decompo- 
sition Z: then any 1-form lu' close enough to lu induces on Ml a rank-2 Zorich 
decomposition Z 1 and this decomposition is homotopic to Z. In particular, all 
such decompositions share the same soul, i.e. the soul is a locally constant 
function of the pair (Ml , lu) . 

Proof. The leaves at the boundary between Mi's genus-1 components of periodic 
leaves Ni and the cilinders of closed leaves are themselves closed and therefore 
stable under small perturbations of w's direction; consequentely, no cilinder of 
closed orbits will disappear in a whole neighborhood of lu £ MP 2 . Since no two 
leaves of a foliation can intersect, this means that open leaves are bounded to 
genus-1 components of A4, homotopic to the Z ones, even for all 1-forms close 
enough to a;. 

Finally, since the homology class I of these rank-2 genus-1 components is 
integer and they change continously, this I must be the same (modulo sign) for 
all them. □ 

The soul I of a Zorich decomposition is a fundamental invariant since its 
knowledge is enough to describe the asymptotic behaviour of the open leaves. 
Indeed, the fact that [Mi] = ±Z means that, in R 3 , the lift Mi of Mi lies between 
a pair of parallel planes perpendicular to I (seen as a direction in R 3 ), so that 
the lift to R 3 of an open leaf, namely an open intersection of Mi with a plane 
perpendicular to ui (seen as a vector in R 3 ), is a curve contained in a plane strip 
of finite width; Dynnikov |Dyn92| showed how these conditions are enough to 
conclude that the leaf is actually a finite deformation of a straight line whose 
direction is the axis of the strip, namely a u> x I" . 
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Since the foliation induced by u> is determined just by its direction H and the 
soul itself can be interpreted as a direction in R 3 , once we fix a surface M we can 
think the soul application as a locally constant application soukx : MP 2 — ► RP 2 . 
We will show in next section that, for a generic polyhedron M, this map is well 
defined on the whole projective plane except for a set of meaure zero ("ergodic 
directions") and its image amounts to a finite number of points. 

Definition 9. The non-empty level sets T>i(A4) = soul^|(Z) corresponding to 
non-zero values of I are called "islands" or "stability zones" of A4. The union 
of all stability zones S(M) = UiT>i(A4) is called the "Stereographic Map" (SM) 
ofM. 

Strange as it seems, this whole construction is the natural model for a very 
concrete physical phenomenon, namely the magnetoresistance in normal metals, 
with the periodic surface being the "Fermi Surface" of a metal, the 1-form 
a strong constant homogenous external magnetic field and the Fermi Surface 
leaves the orbits of the momenta of the metal's quasi-electrons. The presence of 
open leaves is detectable experimentally, and so is each stability zone (provided 
at least two points of it are measured NM98 ) . Experimental plots of the SM 
have been produced in Sixties and Seventies for about thirty metals but only 
recently the first two SM, relative to the Fermi Surfaces of Au and Ag, have 
been reproduced theoretically from first first principles Leo04b] ILeo05bj . 

2.4 The generic (3-irrational) case 

In the previous section we showed that the structure of the foliation in a neigh- 
borhood of rational directions is rather simple: either all leaves are closed or 
their lift is strongly asymptotic to a straight line |Dyn97| , i.e. it is contained in 
a finite-width plane strip and crosses it from one end to the other. Still, this is 
far from being enough, since the sole density of an open set does not preclude 
its measure from being small. 

We will prove in this section that the one above is nevertheless the generic 
situation, namely: 

Theorem 4. The set of directions in RP 2 inducing Dynnikov decompositions 
of genus bigger than one on a generic polyhedron M. has measure zero. 

To prove this theorem, Dynnikov studied the structure of the foliation in- 
duced by a 3-irrational 1-form to on 1-parameter family of surfaces M e — / _1 (e), 
where / is a triply periodic Morse function such that for almost all values of / 
no more than one critical point of to lies on the same leaf. Here, we will repeat 
Dynnikov's steps assuming / to be a generic triply periodic PL function, so 
that almost all of its level surfaces are embedded polyhedra, satisfying the same 
genericity condition with respect to lo. 

The following two lemmas by Dynnikov extend with no change to the poly- 
hedra case: 
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Lemma 2. Given a triply periodic PL function, the set of values for which 
a constant 1-form lu induces open leaves filling rank-2 components is a closed 
connected interval [ei(u), e2(w)]. The functions e and E are continous on the 
whole MP 2 . 

Definition 10. The Dynnikov index w of a critical point c of lu = h a dp a with 
respect to an oriented polyhedron Ai is the product of the "Hamiltonian" index 
of the critical point (+1 for centers and — 1 for saddles) times the sign of the 
scalar product between H = (h a ) and any of the normals to the faces adjacent 
to the critical vertex. 

Lemma 3. The curve r ) u) j(e) — ^ i WiCi(e) is a well defined loop in T 3 and the 
quantity 

J — oo 

is equal to the density of closed leaves on any leaf induced by lu onT 3 . 

Moreover, if {hj} (resp. {h^}) are the heigths of the "positive" (resp. "neg- 
ative") cylinders of closed leaves on Ai, namely those that contain points with 
smaller (resp. bigger) values of f , it turns out that it is possible to choose M 3 
representative Ci of the critical points so that 

x u j(e) = J2 <H ' w M e ) >= Y h i ~ Y h k 

i 

In the weighted sum of all critical points are also contained the pairs center- 
saddle, that can be removed by homotopy and have nothing to do with the 
topology of the foliation. If we do not include them in the summation, we are 
left with a new quantity that tells us the density of non-trivial (in M.) closed 
leaves and therefore is able to spot whether there are cylinders of non-zero 
heights. In particular, to have a full ergodic situation, i.e. a leaf dense on the 
whole M., this function must necessarily be zero. 

Lemma 4. Be {c* op } the subset of the topological saddles of the pair (u>,f). 
Then the reduced curve 7^°y(e) = 'Yli i WiC t ° p {e) is also well-defined in T 3 and 
the reduced Euler density ~X° P f(e) = i* top w is equal to the density of closed 

leaves non homotopic to zero in M. . 

The function ~X° V f is a strictly increasing non continous PL function with 
respect to e; its discontinuity points are exactly the non-proper values of f. 

Proof. The original curve 7 w ,/(e) is well-defined because w i = |Dyn97| ; 
since the Dynnikov indices of a pair saddle-center are opposite, the restricted 
sum J2top w i ^ s s ^ m zer0 an d therefore also 7^°^(e) is well defined. 

Considering only the sum of the topological saddles is equivalent to cancel 
from a leaf all closed leaves that are homotopic to zero in Ai, so that the averaged 
Euler characteristic now is only relative to the non-trivial closed leaves only. 

Now consider a positive cylinder: since the values inside are smaller with 
respect to the values on the cylinder, the height of the cylinder increases with 
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e; for the same reason, negative cylinders decrease their height. Since xj'j(e) = 

Stop ~ Stop ^fc ' then the function xL°/( e ) ^ s strictly increasing in its conti- 
nuity points. The function fails to be continous when new pairs of topological 
saddles are created or destroyed, namely in the non-proper values of /, and it 
jumps exactly by the height of the cylinder(s) created or destroyed. □ 

Corollary 2. If lo induces on A4 e a Dynnikov decomposition of genus bigger 
than 1, then c\{lo) = e^iuS) = eo, i.e. lo induces only closed leaves at any other 
level. 

Corollary 3. At almost all levels of f the measure of the "ergodic" directions 
is zero. 

We have now all ingredients to prove theorem 0] 

Proof. We will just discuss the case of "full ergodicity" , namely the case of 
directions giving rise to leaves dense on the whole polyhedron; the same line of 
arguments extends to the lesser ergodic cases. 

If lo induces on M. fully ergodic leaves, then of course no cylinder can appear 
and therefore, for any Morse PL function / such that M = / _1 (0), it must 
happen that x^°y-(0) = 0. Let us consider now x* op as a function of lo and 

e: then the surface X = (x^° p ) _1 (0) C MP 2 X R, for a generic function /, is 
transversal to the sections MP 2 x {e}: indeed in a projective chart, say h z = 1, wc 
have that d^x^ P f = Si w i( c l° P ) x an d similarly for h v , so that the points on X 
where the gradient is zero are exactly the points where 7^°^(e) = J2i w i c i° P = 0- 
This condition is non-generic, that finally proves the claim of the theorem. □ 

2.5 Structure of the Stereographic Map of surfaces and 
triply periodic functions 

The results of the previous two sections show that the SM S(A4) of a generic 
surface is the disjoint union of a countable set of open sets Vi(M) ("islands"), 
each labeled by a I 6 PH2CT 3 , Z), immersed in a sea of directions that give rise 
only to closed leaves. According to our intuition of the system and the numerical 
experiments made to date, we conjecture that generically the number of islands 
is finite, but no rigorous proof of this fact exists. As matter of fact, this structure 
is exactly the one guessed, from symmetry consideration, by the physicist I.M. 
Lifschitz and his Karkov school about fifty years ago LP59l lL"P60| . 

The boundaries of the islands are reached when the last pair of genus- 1 
components collide because of a cylinder collapse and therefore are charaterized 
by the presence of (at least) a pair of inequivalent critical points on the same leaf. 
The set of these directions is the countable union of the curves < H, Ci~cj >= 0, 
i 7^ j; such curves in the polyhedra case are all straight lines, so that every island 
is actually a (not necessarily convex) polygon. There are reasons to believe that 
these polygons are convex for low genus, i.e. at least for genus 3 and 4, but it 
is easy to build examples of high-genus polyhedra with islands that are either 
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non connected or connected but non convex or even connected but non simply 
connected. 

A crucial observation by Dynnikov |Dyn97| allows to associate a SM also to 
triply periodic functions: 

Theorem 5. Be f a Morse triply periodic function: then, if lu induces on 
M.eO a Zorich decomposition Z e Q, it induces a Zorich decomposition Z e for all 
e S (ei(w), e2(w)) and all these decomposition share the same soul. 

Definition 11. The island T>i(f) corresponding to the label I is the union of 
the corresponding islands of f 's level sets: T>i(f) = U ee sT> i(A4 e ). The SM of f 
is the union of all its islands: S(f) = U;2?;(/) 

The SM corresponding to functions are generically drammatically different 
from those corresponding to surfaces. 

First of all, since rational directions induce necessarily Zorich decomposi- 
tions, the set of islands S(f) is now always dense in RP 2 , and of course lemma[21 
also tells us that in this case the sea of directions giving rise to closed leaves only 
dried up, since every direction either belongs to an islands (or its boundary) or 
is ergodic. 

Moreover, the following property shows that the islands can be sorted in a 
rather complex way: 

Theorem 6. Generically every two zones meet transversally and in a countable 
number of points. 

Proof. No point belonging at the same time to two different zones T>i t and T>i 2 
can have irrationality degree bigger than 2, since it must contain the integer 
direction l\ x l 2 . In the smooth case this would be enough, since the boundaries 
are smooth curves and they generically contain only a countable number of 2- 
irrational points and no rational point. In the PL case though the boundaries 
are actually segments of straight lines and threfore they are actually contained 
in the set of the directions with irrationality degree smaller than one. 

Nevertheless, the theorem holds also in this case for the following reason: 
since every direction at the boundary between two zones is perpendicular to the 
direction l\ x l 2 , then their set is the straight line (in MP 2 ) passing through l\ 
and h- Generically none of the two labels falls on the boundary and therefore 
two zones can meet in a number of points not bigger than the number of sides 
of the island with the smaller number of sides. □ 

Corollary 4. Either there is a single zone, i.e. it exists a label I such that 
T>i(f) = RP 2 , or there are countably many zones and they are dense in the 
whole projective plane. 

Since the islands meet trasversally, the non trivial SM will look like 2- 
dimensional Cantor sets. Retrieving numerically such fractals is not trivial 
since it involves, in general, to analyze the system at several values of /, but 
there is a class of interesting (non generic) cases in which it is actually enough 
to analyze a single level surface of / to get the entire fractal picture: 
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Theorem 7. Be M a polyedron whose interior is equal to its exterior, modulo 
the group G ~ R 3 x Z2 of translations and inversion of the three axes. Then 
S(A4) — S(f) for any function having M. as (connected component of a) level 
set. 

Proof. By symmetry, we can build a function / such that M. = / _1 (0) and that 
/ _1 (e) is equal to e), modulo G. Since also the bundles of parallel planes 

are invariant by G, it turns out that for such an / the interval I(u>) of existence 
of open orbits relative to any 1-form lu is of the form I{to) = [— a(u>), a(co)] 
and therefore G I(uj), G MP 2 , i.e. at the zero level every lj induces open 
leaves. □ 

In particular, all triply periodic functions / such that /(c — x, c — y,c — 
z) = —f(x,y,z) belong to this class, and indeed the only fractals analyzed 
numerically to date are relative to this kind of functions. 

Finally, we cite an important property that ties the set of all labels relative 
to the islands of the SM of a function with the set of ergodic direction [Leo04a 
ILeo05a|: 

Theorem 8. The closure of the set of all labels is the disjoint union of the set 
of all zones boundaries and the set of ergodic directions. 



3 A concrete case study 

As pointed out in the previous section, to date a picture of the fractal has been 
numerically produced for only two functions, an analytical one and a piece- wise 
quadratic one |Leo04aj . 

The analytical one is f(x, y, z) — cos(27ra) + cos(27ry) 4- cos(27rz), invariant 
with respect to translations by integers r = I? cl 3 , that gives rise to genus- 
3 level surfaces in the range (—1,1) and spheres at every other non-critical 
level. This function represents the simplest non-trivial case possible from the 
topological point of view, since any triply periodic connected surface of genus 
smaller than 3 lies between two parallel planes and therefore the aymptotics of 
plane sections is easily found. Its zero level is rather special: it is known as 
the Schwarz primitive function (or plumber's nightmare) and was studied by 
Schwarz in 1890 as one of the first examples of triply periodic minimal surface. 

From the computational point of view, SV = f _1 (0) has three important 
properties: 

SP1 its interior is a translate of its exterior, so that S(SV) = S(f); 

SP2 it is invariant with respect to the natural action of the tetrahedral group T<j 
on the unitary cube, so that the whole SM can be obtained, for example, 
extending by symmetry to the whole KP 2 the data obtained for the triangle 
with vertices [(0,0,1)], [(1,0,1)] and [(1,1,1)]; 
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Figure 3: The simplest PL approximation of the cosine function. 



SP3 the two cylinders, one negative and one positive, have the same height, 
so that it is enough to examine just one of the four topological critical 
points at the base of the two cylinders to retrieve all information about 
the structure of the foliation. 

The piecewise quadratic function is g(x,y,z) — J2q{x 1 ), where fj is the 
simplest piecewise quadratic function having the same symmetries of the cosine 
function. Its level sets have the same behaviour as the function above but the 
expression of the critical points as function of the direction of the 1-form and 
the level of the function is so simple that allows a comparison between analytical 
and numerical data also at levels different from zero. 

Nevertheless, in both cases the number of triangles needed to describe in 
sufficient detail the surface is so big (between 10 5 and 10 6 ) to make impossible 
to improve the resolution of the results obtained in |Leo04a] . at least until a 
new algorithm is found or some "ad hoc" trick used. 

3.1 The polyhedron 

A natural way to improve the resolution of the numerical analysis of the problem 
is to consider PL functions, since in this case the number of triangles needed to 
describe their level surfaces can be as low as of the order of 10 1 . Moreover, the 
description of the surfaces is in this case exact rather than an approximation. 

The simplest case to study, in order to take advantage of the extremely 
convenient properties evidenced in the SV case, is the PL function f)(x,y, z) = 
where rj is the function shown in Fig. El The polyhedron Vo = f) _1 (0) 
is a PL embedding of a genus-3 surface in T 3 having 8 hexagonal faces, 20 edges 
and 12 vertices; the smallest triangulation for Vq takes 4 • 8 = 32 triangles, 
32 • 3/2 = 48 edges and 12 vertices, that gives the expected Euler characteristic 
x(Pa) = F - E + V = -4 = 2- 2-3. The basic cell of the lift V Q c M 3 
in the unit cube is a truncated octahedron (Fig. ; it is noteworthy to notice 
that, like its smooth analog, also this surface is, in the discrete sense, a minimal 
surface Way04| . 
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Figure 4: Plot of the Truncated Octahedron inside T 3 (left) and of part of its 
image in the universal covering. 

Since exactly four edges meet at every vertex, given any 1-form lu in general 
position with respect to Vo only saddles with index — 1 may arise and therefore 
there are always four vertices that are critical for such to. If the point p\ is such 
a vertex, then the remaining three critical points are P2.3 = (1/2, 1/2, 1/2) ±pi 
and P4 = (1, 1, 1) — p\. In particular in our numerical study we sample the set 
of 1-forms w = {h x ' ,W> ', h z ) such that h x /h z € [0, 1] and hy/h z € [0, 1], for which 
the four critical points are (0, .5, .75), (.5, 1, 1.25), (.5,0, .25) and (1, .5, .25). 

3.2 The algorithm 

In order to generate an approximate picture of the fractal, it is enough to pro- 
duce an algorithm able to evaluate the label, if any, associated to a given 1-form. 
Since obviously no calculator can deal with irrational numbers, the numerical 
study will be limited to rational 1-forms; luckily this is not a big restriction, 
since anyway rational directions are dense in every stability zone. 

Note that the algorithm used for the numerical study of the PL case is a 
simplified version of the more general algorithm we developed to study smooth 
surfaces of genus three Lco04a , since in this case we know a priori the position 
of all critical points and moreover we know their position exactly, so that we do 
not have to correct "by hand" the topology of the critical section. 

The basic idea to retrieve the label, as suggested to me by I. Dynnikov, is 
that the soul I £ PH2 (T 3 , Z) associated to the Zorich decomposition Z induced 
by to £ T>i is in 1-1 correspondance with the rank-2 sublattice of ffi(T 3 ,Z) 
obtained as the image, through the map i» : H\(M.,1j) — > -ffi(T 3 ,Z), of the 
open leaves in M that have zero intersection number with the closed leaves 
populating the cylinders of Z. Indeed, every cycle lying on the interior of a 
Nj component of Z has no intersection with the closed leaves that form the 
cylinders Ci, and since all Nj are homologous to each other (modulo sign) the 
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image of all these cycles in T 3 must have rank-2; on the other side, there is an 
obvious 1-1 correspondance between rank-2 sublattices of -ffi(T 3 ,Z) = T ~ Z 3 
and the 2-tori embedded in T 3 = M 3 /r, since every such 2-torus can be spanned 
by a pair of independent rational directions and viceversa. 

The following algorithm N works for genus-3 polyhedra satisfying properties 
SP1-3, in particular for Vo- 

Input: M. - the polyhedron; u> = (l,m,n) 6 Z 3 - the 1-form; x - a critical 
point of u> with respect to A4; tt u ^ x - the plane perpendicular to (l,m,n) and 
passing through x. 

Output: C1.2 - the two critical loops 2 ; h\ t 2 - the homology classes of c\ t 2 in 
M; Hi^ - the homology classes of ci ;2 in T 3 . 
Algorithm: 

Nl retrieve the intersection between M and tt^^; 

N2 check that there are exactly four critical branches and follow them by 
periodicity otherwise exit; 

N3 if no other critical point is met along the path, so that the four branches 
are arranged in a pair of critical loops, store the two loops in the variables 
ci t 2, otherwise exit; 

N4 evaluate the homology class of Ci ; 2 in T 3 and in M (this is actually done 
while executing N2 to speed up the computations time); 

N5 if the saddle is half-open, i.e. if exactly one among is zero, then 
associate to w the complementary h triple, otherwise exit. 

The main outcome of the algorithm is of course the label associated to w. The 
fact that this label is a triple of integers is very important, since an integer 
evaluated numerically with an error smaller than .5 becomes actually an exact 
measure. 

We implemented this algorithm in a C++ library named NTC 3 built over 
an Open Source C++ library named VTK 4 . The choice of the language comes 
from the fact that VTK provides the basic geometric environment and algo- 
rithms needed by the problem, mainly the capability of generating meshes for 
isosurfaces and evaluating intersections between geometric objects. The inheri- 
tance mechanism of the C++ language allows to use transparently all functions 
of a library, hence we used VTK as a starting point and implemented in NTC 
the routines to deal with periodicity and evaluate the homology classes. 

No serious attempt to evaluate the error on such calculations has been made 
to date since no need for it manifested. As check for the reliablity of the result 
were rather used indirect evidences: 

2 since numerically we can study only the 1-rational case, in T 3 the saddles are always 
wedges of circles, i.e. all critical branches close back to the critical point; according to whether 
these loops are or not homotopic to zero in T 3 , their R 3 lift will be open or close. 

3 http:/ /ntc. sf.net/ 
http:/ /www. vtk.org/ 
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• the agreement of the biggest zones with their analytical boundary (Fig. 03), 
obtained through the independent algorithm A; 

• the symmetry of the final picture with respect to the diagonal (Fig. , 
symmetry that was in no way used in the numerical calculations; 

• the agreement of the fractal picture with the labels plot (Fig. I15|l . 

The exploration of the SM was performed in the square [0, l] 2 of the projec- 
tive chart (Lo x /uj z ,ujy/uj z ) by evaluating the label associated to every direction 
at the vertex of a uniform grid of step r and it was repeated for the values 
r = 10 2 , 10 3 , 10 4 . Samplings with r = 10 2 , 10 3 have been succesfully performed 
also for the previous two functions |Leo04a| but the CPU time needed to reach 
r = 10 4 in that case was way too big. It is because of the rather small number 
of triangles needed to describe Vo that the computation became doable. 

3.3 Numerical results for r = 10 2 

This resolution is the lowest one that allows to have a hint of the structure 
of the fractal. About r sections are needed to follow the critical branches for 
a generic direction (to, n,r), m,n G {1, . . . , r}, that takes a time of .5s on a 
1GHz CPU for the evaluaton of a single label and 10 4 x .5s ~ lh for sampling 
the 10 4 directions of the grid (Fig.[JJ|. 

Even from this rough picture it is rather evident a further symmetry of the 
picture, namely the one with respect to the antidiagonal of the square. This 
symmetry does not come, like the others, from the tetrahedral group but it is 
rather of topological nature. The numerical evidence is that, if a 1-form (to, n, r) 
is labeled by L, then its symmetric (1 — n, 1 — m,r) is labeled by L + (1, 1,0) 
but no proof of this fact is known. 

In order to verify the correctness of the algorithm, we found "by hand" 
the analytical boundaries for the biggest zones and compared them with the 
numerical results (Fig. 03). The following algorithm, aimed at the cases similar 
to the cosine one, is a slight modification of the original algorithm introducted 
and used by Dynnikov in |Dyn96| : 

Algorithm A 

Al fix a 1-form u> = (m, n, r) e Z 3 inside some zone (e.g. extracting it from 
the experimental, guessing it from symmetry arguments or simply by trial 
and error); 

A2 retrieve the critical section of ui passing through one of its critical points 
c and make sure it is half open, otherwise exit; 

A3 evaluate the homology class I of the closed critical leaf C; 

A4 rotate u> around some direction till the cylinder of which C is a base 
collapses and identify the critical point c' that is now connected to c 
through a saddle connection; 
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Figure 5: A close-up of the island 2?(2,4,5) (Vo) = D (2,4,5) (I})- Inside the island 
the pairs of critical points at the base of each cylinder are locally constant. In 
correspondance to each sides, there are three different pairings sorted in open 
subsets, labeled in the picture by roman numerals, separated by straight lines 
segments correspoding to directions oj for which the bases of the positive and 
negative cylinders collide, resulting in a saddle connection between two critical 
points. These three segments meet in the single point (.4, .8), that in this case 
happens to be exactly the direction of the label (this property though is not 
generic). In Figure we show in detail the transition between two different 
pairs within this island. 

A5 the equation < H, c — c' >= contains one of the sides of the island; 
follow it in one direction till four critical points fall over the critical closed 
leaf: this is the point when a sides and a new one start; repeat this step 
till the island boundary close up on themselves. 

Since the genus is three, only two cylinders may appear and they will be of 
opposite sign. The pairs of critical points at the base of cylinders are locally 
constant; in the square under investigation the pairings aiepi,pi for the positive 
cylinder and P2,P3 for the negative one, so that boundaries are always given by 
an equation like < H,pi — p4 + L >= 0. See Figures and H3 for a concrete 
example worked out in detail. 

Finally, a picture of the whole fractal can be obtained through the natural 
free action on RP 2 of the tetrahedral group Tj, whose order is 24 (Fig. E{d)). 

3.4 Numerical results for r = 10 3 

This is the highest resolution reached in I. col) la] . In this case for each generic 
direction are needed about 10 3 sections to follow the critical branches, that 
takes a time of 2.5s on a 1GHz CPU for the evaluaton of a single label and 
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Figure 6: Significant examples of critical sections of Vq for ui £ 2?(2.4,5)(7 : 'o)- 
At the "center" of the island (e) there is a saddle connection between the four 
critical points (starting from the highest and turning clockwise) N = p\, E = 
p 2 + (3, -3, 1), S = pi + {1, -3,2) and W = p 2 + (l,-2, 1). In the subzones /, 77 
and 777 (Fig. [SJ the pairs of critical points at the base of the positive cylinder 
are, respectively, N and p4 + (2, —4, 1) (i), S and P3 + (2, —2, 0) (b,c) and E 
and ps + (—1,-1,0) (d, g). The separating segments correspond respectively 
to saddle connections between the pairs of critical points S and W (a), iV and 
S (f) and N and E (h). 
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Figure 7: Numerical plot of the square [0, l] 2 of the SM S(V ) = 5(f)) in the 
projective chart h z = 1 at a resolution r = 10 2 . The color of the islands goes 
from blue to red as the norm of the label grows. In the picture are displayed 
the 106 islands with at least four points out of the total 1741 islands found. 
The missing points that is possible to see in the interior of some of the islands 
are due to failure of the numerical algorithm N, e.g. because of the presence of 
saddle connections. The running time for the 10 4 steps cycle needed to retrieve 
these data takes about lh on a Pentium 1GHz CPU. 




Figure 8: Analytical boundaries of the biggest islands found in Fig. found 
using the algorithm A. All of their boundaries are straight lines segments and 
the corresponding label has been reported, when possible, inside the zone itself. 




Figure 9: Comparison between analytical (Fig[8J and numerical (FigEJ bound- 
aries. 
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Figure 10: Fractal image obtained by letting the tetrahedral group Td act on 
the square in Fig. [7| and then projecting on the disc through the stereographic 
map. 
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therefore about 5 x 10 2 X lh ~ \month for sampling the 10 6 directions of the 
grid (Fig. [Hi)- 

A time of the order of a month is of course rather long, but thanks to the 
diffusion of the Linux OS, and therefore of the possibility to build big Unix 
clusters for cheap using PCs rather than workstations, this is still not too bad 
since it is easy to lower by a factor 10 the computations time just by dividing 
the cycle over as many PCs. This way the running time goes down to just three 
days, that is a rather acceptable time. 

We point out that the situation is radically different in the smooth case: 
indeed in that case there is another variable to consider, that is the resolution of 
the mesh of the surface, that must be increased together with the grid resolution 
to avoid errors in the topology of the curve. The plane sections giving the 
complete intersection between the surface and a 2-torus with homology class 
(l,m,r), l,m £ [0, r], can be as close as 1/r and therefore, if the mesh is too 
rough, there is the risk that the program will jump on the wrong slice. 

Concrete tests show that a mesh resolution of 30, meaning that the mesh 
is produced by dividing the unit cube in a 30 3 uniform grid, is enough for the 
r = 10 2 case but it must raised to at least 60 for the r = 10 3 case, increasing 
the time for a single label evaluation to 15s, an order of magnitude bigger 
than in the PL case. This brings back the time to 3 months for the execution, 
that is indeed the order of the time spent for the r = 10 3 calculations made 
for |Leo04a] . 

From the picture 1111 it is rather evident the symmetry with respect to the 
anti-diagonal. Apart from this, the picture looks qualitatively very similar to 
the pictures found in the previous two cases at the same resolution. 

3.5 Numerical results for r = 10 4 

Increasing by another order of magnitude the resolution, we increase by an order 
of magnitude the number of sections needed to follow a generic leaf, resulting 
in another factor 5 in the running time for a single evaluation of a 1-form label, 
that is now 10s, so that the total running time on a 1GHz CPU reaches 
5 • 10 2 • 20rf = 10 4 rf ~ ZOyears. 

Such a big running time is rather scary and suggests that there is no hope to 
go up by an order of magnitude in resolution without changing some significant 
algorithm step. Nevertheless, this big time can be once again brought down 
to something reasonable by running the code in 20 PCs and by restricting the 
numerical analysis to the upper triangle of the square [0, .5] x [.5, 1] (that reduces 
computational time by a further factor 8). Thanks to all these expedients, the 
running time goes down by two orders of magnitude, reaching about 3 months, 
that is indeed about the time that took to us to collect the r — 10 4 data. 

Note that for the smooth cases this would not be enough, since we must raise 
the mesh resolution to 10 2 and the time for the single evaluation goes up to 50s; 
this, together with the fact that there's a further factor two due to the lack of 
symmetry with respect to the antidiagonal, brings the total running time to 3 
years, definitely not realisticly affordable. 
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Figure 11: Numerical plot of the square [0, l] 2 of the SM S(Vq) = 5(h) in the 
projective chart h z = 1 at a resolution r = 10 3 . The color of the islands goes 
from green to red as the norm of the label grows. In the picture are displayed 
the 1625 islands with at least four points out of the total 10725 islands found. 
The missing points that is possible to see in the interior of some of the islands 
are due to failure of the numerical algorithm N, e.g. because of the presence of 
saddle connections. The running time for the 10 6 steps cycle needed to retrieve 
these data takes about 1 week on a Pentium 1GHz CPU. 
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Figure 12: Box counting evaluation of the Hausdorff dimension of S(Vo) with 
the r = 10 4 data. N n is the number of squares of side 2™ needed to cover 
the complement of S(Vo)> the angular coefficient of the linear fit provides the 
dimension estimate. The same evaluation made with the r — 10 3 data gives a 
very close result, and also restricting the fit by canceling a few points at the 
extremes does not change considerably the estimate. 

In Figures 1 171301 are shown the numerical results, from which it is clear 
beyond any doubt that the SM has a fractal-like self-repeating structure, even 
though no explicit construction is known for it. 

Numerical evaluations of the Hausdorff dimension d-p of the fractal set, 
namely the complement of S(Vo) in MP 2 , has been performed using the box 
counting technique fFig. H2|) and a sort of "area distribution" technique (Fig.ll3|). 

The first, and more reliable, technique involves partitioning the square [0, l] 2 
in 2 n identical squares and evaluating the number of squares needed to cover 
the fractal. With the r = 10 3 data, n can get as big as log 2 (10 3 ) ~ 10 and a 
linear fit gives an evaluation of d-p ~ 1.93. The r — 10 4 data allow n to go 
up to 13, representing the deeper results on the Hausdorff dimension of such 
fractals to date; figure ITSl shows that the scaling law is linear to a high degree 
of accuracy and a linear fit gives again an estimate of d-p ~ 1.93, making us 
rather confident in the accuracy of this numerical result. 

The second one involves counting the number of zones whose size is between 
b n and 6™ +1 , where b > 1. Computations with sevaral small values, between 
2 and 1, were performed and all of them gives rough estimate of d-p ~ 1.76, 
further than expected from the evaluation given by the box counting. It is not 
clear to us the reason of this disagreement, that did not manifest for the smooth 
cases |Leo03j , but it is not impossible that it may be due simply to the fact that 
this evaluation stabilizes at higher resolutions and therefore this evaluation is 
at the current state of things more or less unreliable (as simple numerical tests 
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Figure 13: "Area distribution" evaluation of the Hausdorff dimension of S(Pq) 
with the r = 10 4 data. N n is the number of islands whose area is between 2~ n 
and The estimate of the Hausdorff dimension is given by the double of 

the angular coefficient of the linear fit and varies by ±0.2 by restricting even by 
little the number of points on which the fit is made. 




Figure 14: Logarithmic plot of the islands' areas versus the norm of the cor- 
responding label. There is a very good agreement between the numerical data 
and the Dynnikov conjecture claiming that A(T>i) < C/||/|| 3 . 
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testify). 

Finally, the r = 10 4 data provide enough detail to test numerically a conjec- 
ture by Dynnikov |Dyn99| and give a graphical evidence of theorem|HJ The Dyn- 
nikov conjecture claims that the area of the islands satisfy a relation A(T>i (A4)) < 
C/||Z|| 3 for some positive real number C depending only on the polyhedron.M ; 
the numerical data suggests that the exponent 3 cannot be improved any further 
(see Fig. ITU) . 

According to theorem |H1 the closure of the set of labels of S(Vq), seen as 
points of RP 2 , is equal to the complement of the interior of the islands. At 
this resolution about 5 • 10 5 different labels are found and their image in MP 2 
fFig. I15|) is one of the best indirect checks of correctness of the library NTC. 

4 Conclusions 

We proved in this paper that the structure of foliations induced on triply periodic 
embedded polyhedra by constant 1-forms is identical to the one induced on 
smooth surfaces, and we believe that the same line of arguments can also be 
used to further generalize the theorems to embedded piecewise smooth surfaces. 
This fact extends to "Morse PL functions" the association of a SM; in the 
most interesting cases such SM have a fractal nature (this condition is true for 
an open subset of all triply periodic PL functions, e.g for all triply periodic 
functions close enough to f)). 

Surfaces satisfying property SP1 (see section |3J) are particularly rich, their 
SM being equal to the SM of any function having them as level set. We exploited 
this fact by studying the case of the triply periodic surface T$ obtained extend- 
ing in the three coordinate directions a truncated octahedron. The simplicity 
of the triangulation of the surface allowed us to improve by an order of magni- 
tude, respect to the results obtained in Lco04a , the resolution of the numerical 
analysis of the fractal and, as a consequence, to perform several numerical tests 
on conjectures and theorems. 
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Figure 15: Plot of the labels of the 494041 islands found in the numerical 
analysis, at the resolution r = 10 4 , of the portion of the SM S(Vq) contained 
in the triangle of vertices (0, 1), (.5, 1) and (.5, .5) in the projective chart h z — 
1. According to theorem |H1 the closure of the set of labels is equal to the 
complement of the interiors of the islands; the striking closeness of the two 
pictures is one of the best indirect tests of the correctness of the implementation 
of the N algorithm in our C++ library NTC. 
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Figure 16: The next 14 pictures show, in full detail, the fractal structure of 
S(Vo) = 5(f)) sampled at the resolution r = 10 4 . In order to minimize the CPU 
time, the numerical analysis was limited to the triangle [(0,1,1)], [(1/2,1,1)] 
[(1/2, 1/2, 1)], divided above in smaller squares and triangles by a uniform grid 
and labeled by the corresponding figure number; the full picture can be retrieved 
by applying the "topological" symmetry characteristic of this surface, namely 
the symmetry with respect to the antidiagonal, and then the 24 symmetries 
coming from the action of the tetrahedral group T^. To retrieve these data we 
used about 20 1GHz Pentium III CPUs for about 3 months. In all pictures, the 
color of the islands goes from green to red as the norm of the label grows. Note 
that the square [.4, .5] x [.9, 1] is fully contained inside the island 15(1,2,2) and 
therefore its picture will not be shown. 
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Figure 17: Detail oiS(V ) in [0, .1] x [.9, 1] 




Figure 18: Detail of S(V ) in [.1, .2] x [.9, 1] 




Figure 19: Detail of S(V ) in [.1, .2] x [.8, .9] 




Figure 20: Detail of S(V ) in [.2, .3] x [.9, 1] 




Figure 21: Detail of S(V ) in [.2, .3] x [.8, .9] 




Figure 22: Detail of S(V ) in [.2, .3] x [.7, .8] 




Figure 23: Detail of S(P ) in [.3, .4] x [.8, .9] 




Figure 24: Detail of S(V ) in [.3, .4] x [.7, .8] 




Figure 25: Detail of S(V ) in [.3, .4] x [.6, .7] 




Figure 26: Detail of S(V ) in [.4, .5] x [.9, 1] 




Figure 27: Detail of S(V ) in [.4, .5] x [.8, .9] 




Figure 28: Detail of S(V ) in [.4, .5] x [.7, .8] 




Figure 29: Detail of S(V ) in [.4, .5] x [.6, .7] 




Figure 30: Detail of S{V Q ) in [.4, .5] x [.5, .6] 



